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For diffusive systems that can be described by fluctuating hydrodynamics and by the 
Macroscopic Fluctuation Theory of Bertini et al, the total current fluctuations display uni- 
versal features when the system is closed and in equilibrium. When the system is taken 
out of equilibrium by a boundary-drive, current fluctuations, at least for a particular family 
of diffusive systems, display the same universal features as in equilibrium. To achieve this 
result, we exploit a mapping between the fluctuations in a boundary-driven nonequilibrium 
system and those in its equilibrium counterpart. Finally, we prove, for two well-studied pro- 
cesses, namely the Simple Symmetric Exclusion Process and the Kipnis-Marchioro-Presutti 
model for heat conduction, that the distribution of the current out of equilibrium can be 
deduced from the distribution in equilibrium. Thus, for these two microscopic models, the 
mapping between the out-of-equilibrium setting and the equilibrium one is exact. 

§1. Why studying current fluctuations? 

In one of his 1905 papers, Einstein 1 ) establishes, for the motion of small spheres 
in suspension, a relationship between their diffusion constant and the fluid's viscosity, 
that is a relation between fluctuations in equilibrium and a response coefficient when 
the system is driven away from equilibrium by an infinitesimal force. This instance of 
a fluctuation-dissipation theorem is a particular case of the Green-Kubo relations 2 ) 
which state that, if Q(t) = J j(r,t')dt'dr denotes the total current associated to 
some locally conserved observable, then the variance of Q in equilibrium teaches 
us directly about the transport properties of that very observable. For example, 
if Q is the total particle current in a fluid with mean density p, then its variance 
a = (Q 2 ) c /t verifies 3 ^ 

a = 2Dp 2 k B TK T (1-1) 

where kt = ~§p 1S the isothermal compressibility and D is the diffusion constant. 
That current fluctuations teach us about (slightly) nonequilibrium physics is an 
idea that can pushed forward: nonlinear response coefficients -the so-called Burnett 
coefficients- can be related to higher cumulants of the current. The recent upsurge 
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of interest for current fluctuations can be attributed to the numerical work of Evans, 
Morris and Cohen 4 ) and to the mathematical breakthrough of Gallavotti and Co- 
hen, 5 ) who showed the existence of an extended Einstein's relation applying to the 
entire distribution of the current instead of bearing on its first nontrivial moment. 
This relation takes the form of a particular symmetry property of the distribution 
of Q and now goes by the name of fluctuation theorem. In fact, even for a system in 
equilibrium, current fluctuations tell us about how far the system has been wander- 
ing away from its typical realization. It was indeed recently realized that in order 
to cope with a given current fluctuation, the system may have to adopt a strongly 
heterogeneous configuration. 6 )~ 8 ) 

The goal of the present paper is to identify what the generic properties of the 
current distribution are in systems whose dynamics is diffusive and that can be de- 
scribed by fluctuating hydrodynamics. We shall focus both on systems in equilibrium 
and on systems driven out of equilibrium by boundary constraints. How to exploit 
fluctuating hydrodynamics to obtain predictions regarding current fluctuations (and 
other physically relevant quantities) has been formalized by Bertini et al. into the 
Macroscopic Fluctuation Theory (see also 9 )' 10 ) in another context). The latter, al- 
beit formulated in a physicist's language, will be at the basis of part of the results 
presented here. But the newest and strongest results of this work are concerned 
with a puzzling correspondence between the current distribution in and out of equi- 
librium in two well-studied microscopic models, the Simple Symmetric Exclusion 
Process (SSEP, a model for particle transport) and the Kipnis-Marchioro-Presutti 
(KMP) model for heat conduction. 

Let us now formulate the results we have obtained. Let Q{t) be the total (space 
and time integrated) current flowing through the system over a given time window 
[0, t\. Denoting by j = our interest goes to the distribution of j, which decays 
exponentially with the extensive variable t, as t goes to infinity, 

Probj^ = j\ ~e w W (1-2) 

Alternatively, we shall focus on the generating function of Q, 

Z(s,t) = (e~ sQ(t) ) (1-3) 

which also plays the role of a dynamical partition function for space-time realizations 
of the process in which the current is constrained to adopt a given mean value fixed by 
the conjugate variable s. In this language, the generating function of the cumulants 
of the current, ip(s), given by 

ln(e~ sQ ) 

ip(s) = lim = maxj{7r(j) - sj} (1-4) 

i-^oo t 

can also be viewed as a dynamical free energy (the Legendre transform of the entropy- 
like function 7r(j)) whose physical content is related to the nature of the various 
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dynamical phases able to convey a given current. As for its equilibrium thermody- 
namics counterpart, ip(s) is the quantity to investigate if one wishes to bring forth 
universal features. And indeed, we shall prove that ifi takes a universal scaling form 
for systems in equilibrium. We shall further demonstrate that the same universal 
form holds for boundary-driven systems, with however some important restrictions 
on the phenomenological coefficients D and a. Finally, for both the SSEP and the 
KMP process we will show an exact equality between tjj(s) calculated for the equi- 
librium system and that calculated out of equilibrium. 

This paper is organized as follows. We begin in section [2] by giving two exam- 
ples of microscopic systems that can be described by fluctuating hydrodynamics at 
a coarse-grained scale. Then we explain in section [3] how determining current large 
deviations amounts to evaluating the saddle point contribution of a path-integral. 
This technical step is then put to work in equilibrium (section 0]) and out of equi- 
librium (section E]) where we recall for completeness some of our previous results. n ) 
The new results of this work are presented in sections [6] and [7J there we explain 
how the distribution of the current in the boundary driven SSEP or KMP can be de- 
duced from its expression in the absence of a drive, at the microscopic level. Physical 
conclusions and yet open questions are gathered in section 

§2. Fluctuating hydrodynamics, two examples 

2.1. The Simple Symmetric Exclusion Process (SSEP) 

The Simple Symmetric Exclusion Process can be viewed as model for the trans- 
port of particles on a one-dimensional lattice in which each site can be occupied, at 
most, by one particle. Each particle hops randomly (with a unit rate) to either of its 
two nearest neighbors. The mutual exclusion constraint is the source of all interac- 
tions between particles. Denoting by rii(t') (0 < t 1 < t) the local occupation number 
at site i (a binary variable), we construct an occupation field p(x,r) = n,(i') which 
is assumed to possess smooth variations at the space and time scales x = i/L and 
r = t/L 2 . We refer to 12 ) for an explicit construction of the required coarse-graining. 
It can be shown 13 )' 14 ) that the evolution of p(x, r) is given by the following Langevin 
equation, 

d rP = -d x j, j = -Dd xP + i (2-1) 

where the Gaussian noise £ has correlations (S,(x, t)£(x' , r')) = giifelii S(x — x')5(t — 
t'). The functions D{p) and o~(p) that appear in the expression of the local particle 
current j are given by D = 1 and a(p) = 2p(l — p). 

2.2. The Kipnis-Marchioro-Presutti (KMP) model 

We adopt the formulation of Giardina et aZ. 15 )' 16 ) that describes the Kipnis- 
Marchioro-Presutti model for heat conduction 17 ) in terms of a Langevin process. A 
collection of L harmonic oscillators on a one-dimensional chain are subjected to the 
instantaneous thermal noise produced by their nearest neighbors. Let Xj denote the 
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position of oscillator j, whose evolution is given by 

dx 



3 



dt 



-x 



+ Xj+xrjjj+i - Xj-irjj-ij (2-2) 



where the ltd convention is used and where the are Gaussian white noises with 

variance unity. The coupling of Xj to its nearest neighbors arises through the local 
and fluctuating temperatures x a_ x and x? +1 imposed by its two nearest neighbors. 

In this model, there is local conservation of the energy ej = -4-. Assuming that the 
local energy field has smooth variations at the scales given by x = j/L and r = t' jL? 
(with < x < 1 and < r < t jL? , where t is the macroscopic observation time) , the 
theory of fluctuating hydrodynamics allows us to write that the local energy field 
p(x,r) = Sj{t') evolves according to 

d T p = -d x j, j = -Dd x p + £ (2-3) 

where the Gaussian noise £ has variance t)£(x', t')) = ^felli 5{x — x')5{t — r'). 
For the KMP process of (|2-2p . the functions D(p) and a(p) are given by 

D{p) = 1, a{p) = 4p 2 (2-4) 

We shall not prove this result here and we refer the reader to Bertini et a/. 18 ^ 23 ^ 
and references therein for an introduction to the macroscopic fluctuation theory, and 
to 12 ) for a physicist's approach. 

2.3. General framework 

We now summarize the hypotheses at the basis of fluctuating hydrodynamics. 
The relevant degrees of freedom, be they discrete (as in the SSEP) or continuous (as 
in KMP) are described at a coarse-grained level by a density field p(x,r), in space 
units x = i/L where the system size is unity and the running time is scaled by the 
typical diffusion time at the scale of the system's size, r = t'/L 2 (0 < t' < t). At the 
scale given by the system size, fluctuations are asymptotically small, which accounts 
for the noise in the Langevin evolution equation 



d T p = -d x j, j = -Dd x p + £ (2-5) 
having a variance with a 1/L dependence, 

(£(x, t)Z(x>, t')) = a{p{X L T)) 8(x - x')S(r - r') (2-6) 

The weakness of the noise in the large system size limit is the key ingredient that 
makes our calculations possible, as we shall now present. 

§3. A saddle point calculation 

We start from the Langevin equation (|2-5j) for the field p(x, r) and from the ex- 
pression of the total time and space integrated current Q(t) = I? dr Jq dx j(x, r) , 
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whose generating function we write in the form of a path integral based on the 
Janssen-De Dominicis 24 )' 25 ) mapping: 



Z(s,t) = (e- S ®) 



VpVpe 



-LS[p,p] 



where the action is expressed as 



S 



t/L 2 



dr / dx 



pd T p + D(p)d x pd x p - -a{p)(d x p - sL) 2 - (sL)Dd x p 



(3-1) 



(3-2) 



We denote by p(x, r) = p(x, r) — sLx. As was pointed earlier 12 )' 26 ) the path integral 
in (|3-ip calls for a saddle point evaluation in the large system size limit L —* oo. We 
denote by p c (x,r) and p c (x,T) the solutions to 



d T p - d x (Dd x p) + d x (ad x p) = 



rri / 

-— = d T p + d x (Dd x p) + Ud x pf = o 

op 2, 



(3-3) 



Equations (j3-3[) must be complemented with the appropriate boundary conditions. 12 ) 
To leading order in L the partition function reads 

-LS[pc,p c ] 



Z(s,t) 



(3-4) 



We shall assume that the saddle point solution (p c , p c ) is stationary (this issue was 
discussed e.g. in 6 )' 21 )' 22 )). This assumption, when not fulfilled, is signalled by insta- 
bilities that are interpreted as phase transitions. 6 ) -8 )' 21 )' 22 )' 27 ) Therefore, to leading 
order in the system size we have that 

-l 



saddle 



L 



p(sL) 



dx 



D{p c )d x p c d x p c - ^a(p c )(d x p c ) 2 



(3-5) 



Of course, as in any saddle point calculation, it is important to evaluate the 



fluct 



to the asymptotic behavior given in (!3-5p . This is done 



leading corrections i(^(s) 

by expanding the action 5 around the saddle to quadratic order in the deviation 
from the saddle (j) = p — p c and (ft = p — p c , 



d T (ft + D(p c )d x (ftd x (ft + D'(p c )d x p c (ftd x (ft + ^D"(p c )d x p c d x p c (ft 2 

+ D'(p c )d x p c d x (ft(ft - ^o-(p c )(d x (ft) 2 - a'(p c )d x p c (j)d x (ft - ^a"(p c )(d x p c )(ft 2 

(3-6) 

and by integrating out the resulting quadratic form. Note that the latter step, 
which requires diagonalizing the quadratic form (13-6|) . may prove difficult when the 
coefficients of the quadratic form are space dependent, or, equivalently, if the saddle 
point solution is not homogeneous. In the next two sections, we implement the 
program we have just sketched in two distinct settings: for a closed equilibrium 
system and for an open system driven out of equilibrium by boundary constraints. 
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§4. In equilibrium: closed systems with periodic boundary conditions 



We first consider closed systems with periodic boundary conditions. 8 ) The solu- 
tion to the saddle point equations f)3-3[) is indeed rather simple to find, namely 



p c (x) = p, p c (x) 



-sLx 



(4-1) 



where p (with no argument) is the space averaged density. This leads to p(X) 
^a(p)X 2 , which, with A = sL, also reads ip(s) 



saddle 



L^crs 2 . Corrections to the 



saddle arising from integrating out the quadratic fluctuations around the optimal 
profile p c , p c are not hard to evaluate, since the quadratic form (|3-6p has constant 
coefficients. To do so we expand 4> and <p in Fourier modes indexed with wave vectors 
q = 27rn, with n £ Z, as imposed by the periodic boundary conditions. We find that 
the contribution of the determinant reads 



fluct 2L 2 



-y 

T.2 Z-> 



1 L 



which we rewrite in the form 



Dq 2 



D<r{Dq 2 -^(sL) 2 



^ (S) " ~2T S ~ I? 9 \ iW (sL) 



(4-2) 



(4-3) 



where is a universal scaling function, 8 ) a representation of which is given in terms 
of the Bernoulli numbers -E?2n : 



B2k- 



k>2 



r(k)r(k + i 



■{-2x) k 



(4-4) 



The scaling function & has a branch cut along the positive real axis when x > it 2 /2. 
If the argument x = (sL) 2 of hits the value ir 2 /2 upon varying s this signals 
that the basic hypotheses underlying the saddle point calculation are not fulfilled, 
e.g. that the stationary saddle point solution becomes unstable. 21 )' 22 ^ We refer the 
reader to Bodineau and Derrida 6 )' 7 ) for an interpretation in terms of dynamic phase 
transitions. 

At fixed value of s and in the large system size limit L — > oo, the limiting 
behavior of if) is given by 



ij}{a) 



V2 



1 2 
— OS + 

2 3vr 



CT 3 / 2 | S | 3 



+ o{\s\ 



(4-5) 



whose fourth derivative is singular at s = 0. This was interpreted by Lebowitz and 
Spohn, 28 ) in the particular case of the SSEP, in terms of the Burnett coefficients 
being infinite. This result is shown to apply irrespective of the explicit expression of 
D and a. 
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§5. Out of equilibrium: open boundary-driven systems 

We now turn to an open system with the same bulk dynamics as that given by 
(|2-5p . in contact at its boundaries with reservoirs that impose prescribed values for 
the field: p(0,r) = po and p(l,r) = p\. The saddle point equations (|3-3p must now 
be solved bearing in mind these new boundary conditions. A stationary solution does 
exist, although it is now strongly space dependent. This should not be a surprise 
given that already at s = 0, the optimal profile has a nonzero gradient allowing to 
bridge po to p\. In general, the explicit form of p c {x) and p c (x) is difficult to obtain. 
The function p that appears in the rhs of (|3-5p . as calculated from plugging the 
solution (|3-3p using the new boundary conditions into (j3-5|) is exactly the one that 
Bodineau and Derrida 29 ) initially found in their paper on the additivity principle. 
When D(p) is a constant and cr(p) is a quadratic function of p then the analytics 
somewhat simplify and it can be seen by direct calculation 11 ) that, for D = 1 and 
a(p) = c\p + C2P 2 , the saddle point contribution is given by 

ma) = { ;f J arcsinh ^ 2 2 f ru;> ° (5-i) 

l +^(arcsm v /= ^) for lj < K ' 

where uj{\, po,pi) is the auxiliary variable given by 

w(A,po,pi) = |(1 - e clA / 2 ) (d(p a - e-* x ' 2 po) - c 2 (e" clA / 2 _ 1)popi ) (5 . 2 ) 

For the SSEP, o~(p) = 2p(l — p) and one recovers the known 30 )' 31 ) result (the notation 
z = e _A is used in the formula (2.14) of 30 )), namely 

lo(X, po, pi) = (1 - e A )(e~ A p ~Pi~ (e" A - l)poPi) (5-3) 

For the KMP chain of coupled harmonic oscillators, the variable u is now given by 

w(A, po, pi) = \{2{po - pi) - 4A/J0P1) (5-4) 

The difficulty, at this stage, remains to diagonalize the quadratic form (|3-6|) given 
that its coefficients are space-dependent constants. The eigenmodes are not the 
standard plane waves anymore given that translation invariance does not hold. We 
have not been able to carry out this task in general, but we have found a way to 
bypass this technical step when D is constant and a is a quadratic function of p. By 
introducing two auxiliary fields ■tjj and if) defined by 

4> = {dxPc)' 1 ^ + d x p c i>, 4> = d x p c ip (5-5) 

which we substitute into (|3-6|) . and after extensively using (|3-3|) . we arrive at the 
following expression for S 

S = + J dxdr (i)d T iP + Dd x $d x ij> - p{sL)(dJf - ^V 2 ) (5-6) 

The local rotation of the fluctuation fields (|5-5p has allowed to disentangle the space 
dependence and to find a set of variables in which translation invariance is recovered. 
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The action (|5-6p exactly describes the quadratic fluctuations around the saddle in an 
open system in equilibrium, in which the parameter conjugate to the current is now 
s' = i^jp-. We diagonalize (|5-6p with the help of the Fourier modes {sinqx} q=nn , 
n £ N* consistent with the field being fixed at the x = and x = 1 boundaries. The 
conclusion of this section is that for systems having a constant D and a quadratic a, 
we can actually determine the finite size corrections to the large deviation function 
and we find that 



This is the very same function & that appears here for a boundary-driven open 
system as the one that was found when studying its closed equilibrium counterpart. 
We thus draw the partial conclusion that at least for a subclass of systems described 
by fluctuating hydrodynamics (those with constant D and quadratic a), the current 
distribution displays universal features, and these are the same as the ones observed 
in equilibrium. To reach this conclusion, we have resorted to a local mapping of the 
out-of-equilibrium system's fluctuations onto those of a corresponding equilibrium 
system. 

§6. Exact mapping for the driven SSEP onto an equilibrium system 

Let us consider the evolution operator of the SSEP on a one-dimensional lattice 
with L sites with injection rate at the left (resp. right) boundary a (resp. 5) and 
annihilation rate at the left (resp. right) boundary 7 (resp. 0). The hopping rate 
is set to 1. In the present section, and in the next, we find it more convenient to 
study the statistics of the total current flowing between the final site L and the right 
reservoir, for which we denote the conjugate variable A. It was shown explicitly in 11 ) 
that the formal replacement of A with s(L + 1) in the large deviation function allowed 
to pass from the current from the last site to the current flowing through the whole 
system. Thus we consider the evolution operator of the SSEP with the constraint 
that it has to carry a prescribed mean particle current (enforced by the Lagrange 
multiplier A) between site L and the rightmost reservoir. This evolution operator 
can be expressed in terms of the Pauli matrices crj, crj and cr|, and the raising and 
lowering operators = g(cr? ± «<r|), whose algebra is given by 




(5-7) 



[a z ,a ± ] = ±2a ± , [a-,a + ] 



(6-1) 



It reads 



W L (A) 




(Ti 

2 1 2 




(6-2) 



^(eV L+ + i.£-i) +/5 (e-V L --Ia£-I) 



The parameter A is conjugate to the time- integrated current flowing from site L to 
the right particle reservoir. Let us now consider a rotation of the spins 5j = Rsj, 
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where we write the SO(3) matrix R with a Cayley representation indexed by three 
parameters x, y and z, namely R = (J + A)(I — 4)~ x with 



so that, explicitly. 

1 



4 = \ iz -ix (6-3) 




-x 2 — y 2 + z +1 2i(xy — z) 2{y — xz) 

R = \ Z3 TZ3 — I 2i ( x y + z ) x 2 + y 2 + z 2 + 1 -2i(x - yz) 

-2(y + xz) 2i{x + yz) x 2 — y 2 — z 2 + 1 

(6-4) 

We carry out the rotation of the spins in the evolution operator Wl(A) which appears 
in (16-2(1 and we search for a rotation R that allows to interpret the resulting operator, 
when expressed in terms of the new variables s, as an evolution operator for a driven 
and open SSEP with modified rates a', 0', 7' and 5' , and with a modified parameter A' 
constraining the particle current flowing out of the system. The resulting constraints 
read y = —z and 

a' = (l-*)°~(* + *h (6 . 5) 
1 — x 

y = {z-x)a + (l + z) 1 
1 — x 



4 - z)[l + xe x + z(l - e A )] 5 - (x + z) [x + e~ A - z(l - e _A )] 

1 — x 2 

(z - x) [x + e A + z(l - e A )] 5 + (1 + z) [l + xe~ A - z(l - e" A )] /3 



(6-7) 



_ \^ ^-y L ^-rc -r^i-i c ;j u ~r V x t ^ ^ ~r c ;j M ^ 



and the effective A' verifies 



-A' 



x + e A + z(e 



1 + xe~ A + z(e~ A - 1) 



(6-9) 



It is convenient to rewrite the above conditions in terms of the original densities 
Po = -^pj an d pi = -g^p and in terms of the auxiliary parameters a = and 
b = j^p- These now read, with obvious definitions of the primed quantities, 

P'o = (1 + X 1 )P °" X " Z (6-10) 
1 — x 

1 \ \nn \x + e x + z(l — e x )]pi —x — z 

p'l = (x + e" A - z 1 - e" A I \ ^ (641 

1 — x z 

0! = a (642) 

b' = b (643) 

At this stage, we have simply mapped our evolution operator describing the driven 
nonequilibrium SSEP with parameters (a, 0, 7, 5, A) onto another driven SSEP with 
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new parameters (a', /?', 7', 8', A'). 

We now go one step further and we ask if there exists a rotation (that is a pair 
of variables x and z) such that the primed process is in equilibrium, that is, such 
that the stationary densities p' and p\ at the left and right reservoir are equal 

Po = P'l (6-14) 
Such a condition can never be fulfilled at A = 0, but at A 7^ 0, a solution for x and 
z always exists. We have thus established that the nonequilibrium open and driven 
SSEP can be mapped onto an equilibrium open SSEP at arbitrary density. 

It is interesting that we can exploit the freedom to choose the equilibrium density 
to which the original nonequilibrium process is mapped: density p' = p'\ = \ indeed 
plays a special role for the SSEP, since, at this very density, whatever the forcing 
strength A', the density profile remains flat at a value 1/2 at the macroscopic level. 
This makes the computation of the current large deviation function, in equilibrium 
at density \ particularly easy. The condition p' = p'\ = \ leads to 

e- y = (v^ + VTTu) 2 , u = (1 - e A )(e~Vo - Pi - (e" A - l)p Pi) (6-15) 
and hence 

ip L (\;p ,pi,a,b) = tp L ^-21n(V^ + \/T+~tJ); ^, ^,a,b\ (6-16) 

We know that in the large system size limit, ipL (A'; ^, ^, a, 6) = which immedi- 
ately allows us to recover the result of, 30 ) 

1 2 

tpL(^,Po,Pi,a,b) = — (arcsinh^/w) (6- 17 ) 

Li 

We have therefore shown that the cumulant generating function of the current out 
of equilibrium can be inferred from that in equilibrium. 

Moreover, an equality analogous to (|6-16p holds for the full operator of evolution, 
which implies that at fixed rates a and b, the partition function Z(s,t) = (e~ s Q) 
depends on po, p\ and A only through the variable to, for all time t and size L, a result 
in the spirit of. 30 )' 32 ) Last, the exact mapping of this section directly translates at the 
level of the hydrodynamic fields p, p of section[3j Indeed, following, 12 ) the action (|3 2j) 
may be recovered from the evolution operator through the correspondence S + = 
(1 — p)e p , S~ = pe~ p , S z = 2p — 1. The rotation corresponds to the change of fields 

, e-P {x + y+(z- l)eP) {(p - l)e^(x - y) - (z + l)p) 
P = 1 TT~2 2 6-18 

e? = x + ^ + (z ~ 1)e " (6-19) 
e'P(x-y)-z-l K ' 

One checks by direct computation it leaves the bulk action invariant, while the 
boundary conditions become p'(0) = p' , p'{^) = p'\ for x,y,z solution of (I6T0|) . 
(|6-10p . Choosing p' = p'\ = \, one checks that this change of fields becomes (|5-5p 
for the fluctuations around saddle. 
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§7. Exact mapping for the driven KMP onto an equilibrium system 



We consider a microscopic version of the KMP process described in subsection 
(|2.2p where the leftmost (resp. rightmost) oscillator is coupled to a heat bath at 
temperature To (resp. Ti) with an exchange rate 70 (resp. 71). The bulk dynamics 
given in (|2.2p is unchanged but the contact with the heat baths is now described by 

dxi 



dt 
dx L 
dt 



xi + x 2 m,2 - V 270T0C0 

^ + 71 J X L + \/ 2 7i r iCl - ^L-\f]L-X,L 



(7-1) 

(7-2) 
(7-3) 



where £0 an d £1 are Gaussian white noises with unit variance, 
flowing from oscillator L to the bath on the right hand side is 



JL+1 



71 Ti - y/2 r y\Tix L £i 



The heat current 



(7-4) 



where the ltd convention is used. The Fokker-Planck evolution operator for the 
KMP process not only contains the contribution given by 12 )' 15 )' 16 ) that describes 
the unconstrained dynamics, but it also contains A-dependent contributions that 
constrain the trajectories to carry a given mean current whose value is tuned by 
that of A. We find that 

L-l 



Wl(A) = ^2(KjJK j+1 + 1/4) + 70 



K{ + 2T K{ + - 



(7-5) 



+ 7i 



where Kj 



K z L (l - 2AXi) + 2T L Kl + 2X(XT 1 

1/2, K'i 



i)k: 



1 

+ 2 



dj( Xj - 



-k: 



x*/2, KJ = d]/2, Kl >; ... ;) 

The A-dependent terms in (|7-5p can be found directly from a Kramers- 
Moyal expansion as deduced from the Langevin equation for (jT- 3[) and from the 
expression of the current (|7-4p . These operators verify the so-called SU(1, 1) algebra 
relations 

[K z , K^] = ±2K±, 
The SO(2, 1) metric matrix J has elements 



[K,K- 



K z 



(7-6) 



J 



(7-7) 



For KMP, we search, as explained in, 33 ) for a Cayley representation of £0(2, 1) 
isometries in the following way. We search for a matrix A verifying A T J + J A = 0. 
Such a matrix takes the general form 



(7-8) 
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so that the matrix R can now be cast in the form R = (I + A) (I — A) 1 , namely 

- / 1 + x 2 +y 2 + z 2 2i(z + yx) 2(y - xz) 

R = 2 2 2 — 2i-z + 2iyx 1 — x 2 — y 2 + z 2 —2izy + 2ix 

1 + x - y - z y 2(y + xz) 2i(zy + x) 1 - x 2 + y 2 - z 2 

(7-9) 

To each matrix R of SO(2, 1) one can associate a <SC/ (1, 1) transformation that leaves 
the K a, s algebra (|7-6p invariant. This allows us to now proceed along the lines of 
the reasoning carried out for the SSEP. We define kj such that Kj = Rkj and we ask 
whether the evolution operator (|7-5p . when expressed in terms of the new operators 
kj, can be interpreted as the evolution operator of an open and driven KMP process 
with modified bath and current-forcing parameters, j'q, Tq, j[, T[ and A'. This is 
indeed the case provided y = — x and 

. Tn(\ — z) — 2x ,_ 

TL = '- 7-10) 

u 1 + z K ' 



r , = - z -t- zxa) - txjjyi - z -t- ±xa) ^ 711 ^ 



(Ti(l - z + 2a? A) -2a ))(l - z + 2xX) 
T^V 2 

70 = 70 (7-12) 

71 = 71 (7T3) 



and the new A' is given by 

A' = {1 + Z) \ (7-14) 
1 - z + 2xX y 1 

Note that the conditions 7g = 70 and 7^ = 71 are analogous to the conditions 
a' = a and b' = b in the SSEP. There always exists a solution for x such that the 
transformed dynamics describes current fluctuations in an equilibrium system, that 
is with Tq = T[. The latter temperature is then parametrized by z. For each value 
of (x, z), the combination uj = A(To — T% — XTqTi) is left invariant by passing to the 
primed variables, but we have not been able to exploit this fact to recover, by simple 
means, the result (|5T|5-4p . Just as was the case for the SSEP, the A — > limit is 
singular and the mapping fails to hold in that limit. One checks however that (for 
instance imposing Tq = T[ = 1), at fixed 70, 71 the spectrum of the operator depends 
on Tq, T\ and A only through the variable u: 



S P W L (A; T , ^,70,71) =SpW L (-21n(v^ + VT+^);l,l,7o,7i) (7-15) 

a result similar to that of the SSEP (section [6]), which seems to endow u with a 
physical meaning yet to uncover. 

§8. Open issues 

It is well-known 34 ) that boundary driven systems develop long-range correlations. 
It is thus, at first sight, rather puzzling that a local mapping such as the one of 
section ([6]) or (|7|) allows to map a nonequilibrium situation onto an equilibrium one. 
When constraining the dynamics to carry a prescribed mean current imposed by a 
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Lagrange multiplier A, in the long time limit, the physical states associated with 
a given value of A do not display long range correlations. This can be seen by 
combining the explicit evaluation of correlation functions, as done by Bodineau et 
al. 27 ^ with the results of Imparato et al., 11 ^ which gives a finite correlation length 
£(s) = jj-(aa")~2 . Since the long-rangedness disappears at nonzero A, it may be 
less surprising that a local transformation does the trick. In the limit A — > 0, the 
correlation length £(s) becomes infinite which restores the long-range correlations of 
the unbiased dynamics. Such a simplification did not occur in 12 )' 35 ) where density 
large deviations were considered in the absence of a A-drive, which may account 
for the nonlocal transformations needed in that work to map the nonequilibrium 
dynamics onto equilibrium dynamics. 

We do not doubt that similar transformations can be found at the level of fluc- 
tuating hydrodynamics (beyond quadratic fluctuations) for systems belonging to the 
same family as the SSEP and KMP (with a constant diffusion constant D(p) and a 
quadratic noise variance o~(p)). It would be interesting to see the explicit form of the 
continuum analog of our (pseudo)rotations. But of course, a much more interesting 
issue is whether our conclusions hold irrespective of the particular form of D(p) and 
a(p). But that's another kettle of fish. 
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